Trellises for stabilizer codes: definition and uses 
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(Dated:) 

Trellises play an important theoretical and practical role for classical codes. Their main utility is to 
devise complexity-efficient error estimation algorithms. Here, we describe trellis representations for 
quantum stabilizer codes. We show that they share the same properties as their classical analogs. 
In particular, for any stabilizer code it is possible to find a minimal trellis representation. Our 
construction is illustrated by two fundamental error estimation algorithms. 



Introduction. Since the discovery of efficient quan- 
tum algorithms for solving hard classical problcms, many 
cfforts have been devoted to building quantum proccssing 
devices. While small scale prototypes are readily avail- 
able, scalability remains a practical issue because of the 
extreme sensibility of quantum devices to external noise. 
Fortunately, theoretical advances, such as the discovery 
of error correction schemes and fault-tolerant implcmcn- 
tations, have notably cleared the future of quantum com- 
puting. One way of building quantum codes is through 
the stabilizer formalism An (n,k) stabilizer code 

protects k qubits by encoding thcm into an n-qubit reg- 
istcr. The error recovery procedure involves the mea- 
surement of a syndrome (i.e., a vector in F2~ fc ) which 
is used to partially discriminate the actual error E from 
all possible ones. The role of the error model, and in 
general of any a priori information about E, is to allow 
further discrimination in order to find the most likely 
guess E. While one can imagine computing the likeli- 
hood of all possible errors compatible with the measured 
syndrome, such method is impractical when codes are 
large. This is because there are of order 2 n+k such el- 
ements. This problem is well known in classical coding 
theory, but was rightfully ignored in the quantum case as 
block codes that can be physically implemcnted have ex- 
tremely small length. However, with the development of 
quantum communication and the advent of other coding 
strategies, this shall no longer be the case. 

In classical coding theory, one often refies on a graph- 
ical representation of the code, called a trellis, to per- 
forin error estimation. For instance, trellises yield many 
complexity-efficient error estimation schemes for memo- 
ryless channels as well as means of estimating the noise 
parameters. In particular, it can be used to calculate 
with linear complexity the most likely error for a convo- 
lutional code of bounded memory over any memoryless 
channcl. In this article, we apply general rcsults con- 
cerning group codes for classical communication [2j to 
show that a similar representation is available for quan- 
tum stabilizer codes. Two error estimation schemes that 
cxploit trellises are introduced for memoryless channels. 
We show that the complexity of these algorithms is re- 
lated to the number of trellis vertices, and provide a con- 



struction of a trellis which minimizcs this quantity. Onc 
of these algorithms achieves the performance of an algo- 
rithm for convolutional codes that was proposed in jj], 
but with a significantly lower complexity. 

Stabilizer codes: elementary facts. In the rest of 
this paper, some familiarity with quantum computation 
is assumed. This section provides a brief introduction 
to stabilizer codes, for a more detailed introduction the 
reader is redirected to 0, 0, and references therein. 

In what follows, without loss of generality, the quan- 
tum register of interest has n physical qubits. 

Preliminaries. Stabilizer codes rely heavily on prop- 
erties of G n , the n-qubit Pauli group. This group is de- 
fined in terms of the Pauli matrices for a single qubit: 

1= (h°i);X = (Vo),y = (° o'M = (à-°i)- The 
group Gn is the multiplicativc group gencrated by the 
n-fold tensor produets of single qubit Pauli matrices. 

For our purpose here, phases are irrelevant and it will 
be more convenient to work with the effective Pauli group 
G„±g n /{±l® n ;±il® n } (see 0). The elements of d 
will be denoted by I =[I], X =[X], Y =[y], and Zé\Z\. 
Here, [P] denotes the equivalence class of V G Gn, that 
is {±V,±iV}. Note that G n is Abelian, so that we will 
use the additive notation for its group operation. Since 
G n = G™, we often view P G G n as an n-tuple (P i )"_ 1 
with entries in G\. 

The crucial fact about Gn, is that any pair of elements 
V, Q either commutes or anti-commutes. This leads to 
the definition of an inner product "*" for elements of 
G n such that {P l ) t * (Q 4 ), = J2i pt * Q l m °d 2. Here, 
P l *Q l = 1 if P i ± Q\ P i ^ I and Q l ^ I; and P l *Q l = 
otherwise. One can then check easily that P,Q G Gn 
commute if and only if [P] * [Q] =0. 

Error model. Stabilizer codes can accommodate for 
a broad class of channels. For simplicity, only memo- 
ryless Pauli channels will be considered, although the 
tools presented in this paper extend to other memory- 
less channels. Memoryless Pauli channels act on the 
whole n-qubit register as u — > ^(c) = (í* 1 ® V P 2 <8> 
. . . (£> *B n )(cr). Above is a 1-qubit channel whose ac- 
tion on p, a single qubit density operator, is given by 

= Y.£e{i,x,y,z} Vï i(l£\)£p£> wit h Pr»(-) a proba- 
bility distribution over G±. 
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Definition of the code subspace. The code subspace C 
of an (n, k) stabilizer code is the largest subspace stabi- 
lized by the action of S, an Abelian subgroup of Q n . For 
the code to protect k qubits using n, i.e. to be of rate k/n, 
S must be generated hy n — k independent operators Sj , 
and be such that —X® n ^ S. The code subspace is equiv- 
alently defined by n — k eigenvalue equations: \ip) € C if 
and only if Vj G {1, .. .,n— k}, Sj\ip) = 

To study the main properties of these codes, phases are 
again irrelevant. More precisely, it is sufRcicnt to rcprc- 
sent the set of generators of the stabilizer group {Sj}j 
by the set of equivalence classes {Sj} where Sj = [Sj] 
which generate a subgroup S of G n . Using a slight abuse 
in terminology, we also call S the stabilizer group of the 
code and {Sj}j the stabilizer set of the code. 

Quantum convolutional codes. Following the defini- 
tion of 0,0, an (n, k) stabilizer code is convolutional 
with paramctcrs (77, k) if there exists a set of generators 
{Sj}j of its stabilizer group S with an 77-qubit shift in- 
variance property. More precisely, the vàlues Sj must 
be equal to the entries í/j G G\ of an infinite matrix H 
which satisfies H % ^l_ K = H* for every i, j. 

For a convolutional code givcn by the stabilizer set 
{Sj}j, the memory is defined as m = maxj J^ . $(Sj). 
Herc, $(Sj) = if eithcr = I for all k > i or S$ = I 
for all k < i; and 'Í(Sj) = 1 otherwise. It is important to 
note that the memory depends on the set {Sj}j and is 
thus not an intrinsic property of the code. 

Error estimation. The goal of error estimation is to 
infer channel errors from their action on the state of the 
quantum register. In the context of stabilizer codes, the 
necessary information is provided by the measurement of 
the Hermitian operators Sj . 

Let £ G Q n be the actual, yet unknown, quantum 
error that affected the state G C. The measure- 
ment of the operators Sj on |-0') = £|'0) defines a bi- 
nary vector of length n — k called the syndrome of £: 

=(«>(£))"=* = è(i - mSiW)) = 

Among all possible error operators f £ Ç„, only some 
of them are compatible with s(£) (i.e., thcy satisfy 
sÇF) = s(£)). Thcrcforc, they belong to a coset of 
N(S), the normalizer of S in Q n . Equivalcntly, these 
compatible errors T are such that \T\ belongs to a coset 
of S 1 - ={P e G n : P * Q = 0, VQ G S}. While, the 
knowledge of the syndrome restriets the class of errors 
that could have happened during the transmission, the 
error model further discriminates between these elements 
by assigning them probabilities. Error recovery then uses 
these probabilities to find a best guess £ for the actual 
error £ . For instance, maximum likelihood error estima- 
tion consists in finding a most likcly £ compatible with 
the measured syndrome s(£). Trellises are both aimcd 
at computing these probabilities and at choosing a best 
guess efficiently. 

Trellises for stabilizer codes. Considcring the pre- 
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FIG. 1: Trellis representation for the 4-qubit code and syn- 
drome (0,0). Here, |V | = |V 4 | = 1 and |Vi| = |V 2 | = |V 3 | = 
4. 

vious remarks about error estimation, it is natural to 
seek a representation of cosets of S 1 - in which it is easy 
to search for their most likely element. This is the main 
motivation for the definition of trellises for quantum as 
well as for classical codes. In a broader context, this mo- 
tivation extends to the calculation of quantitics on which 
error estimation is based, e.g. the likelihood function, the 
a posteriori qubit error probability, etc. 

Definition. In the rest of this section, S denotes the 
stabilizer group of a quantum code with parameters 
(n, fe), and {Sj}"~^ its stabilizer set. 

An n-section trellis rclativc to the stabilizer set 
{Sj}™Zi and syndrome s e is a directed graph with 
the following properties: 

(1) Its vertices can be grouped into n + 1 sets V^, with 
|Vo| = |V„| = 1. The set is called the i-th state space 
of the trellis 

(2) Its edges are directed and can be grouped into n 
sets An edge e=(vw) is said to be issued at vèrtex 
v and ending at vèrtex w. Edges in are issued from a 
vèrtex of and end at a vèrtex of V^. The set Ei is 
called the i-th section of the trellis. 

(3) An edge e e E^ bears a label l(e) such that l(e) € 
G\. We say that l(e) is the Pauli-label of e. 

(4) Each element P € G n with syndrome s is asso- 
ciatcd to a unique directed path (e 1 , . . . ,e n ) such that 
l(e l ) = P\ 

An example of trellis is considered in Figure ^ for 
the 4-qubit code with stabilizer group generated by 
{XXXX, ZZZZ}. 

Construction. Given a code stabilized by S, there are 
many possible trellises representing the coset of S 1 - of 
syndrome s. Wc provide here a simple construction for 
which the number of vertices in each is bounded by 
2 n ~ k . In analogy with classical codes, this trellis will be 
called the Wolff trellis of the code defined by the stabi- 
lizer set {Sj}j relative to the syndrome s. 

For every i G {0, let ttí be a mapping 

from G n to itself defined by tt^P 1 , P 2 , . . . , P") = 
(P 1 , . . . , P\ I) (with the convention that 

TV Q (P x ,P 2 ,...,P n ) = (/,...,/)). Let P s be an ar- 
bitrary, but fixed, element of G n with syndrome s. For 
every i S {0, . . . , n}, is a subset of F™- fc defined by 
{(^•★TT^P))^ : P G P s + S- 1 -}. A vèrtex v G V, is con- 
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ncctcd to vcrtcx w G Vi + i with an edge labelled by E iff 
there exists P E P s + S 1 - such that v = (Sj * 7Ti(P))"~*, 
w = (Sj *7r l+ i(P);=í and P l+1 = E. 

An exemple of trellis obtained in this way is given in 
Figurc[2]for the 5-qubit code associated to the stabilizer 
set {ZXIII,XZXII,IXZXI,IIXZX}. 

As hinted by this examplc, the Wolff trellis of convo- 
lutional codes has state-spaces that are bounded in size 
by 2 m , where m is the memory associated to the chosen 
stabilizer set {Sj}j. Note that this bound is achieved for 
all indices i such that m = J2j tK^])- 

Minimality. As it will be seen below, the complex- 
ity of many useful algorithms using trellises is linear in 
their number of vertices. This raises the issue of finding 
a trellis which minimizes this quantity. If one is will- 
ing to change the stabilizer set for the code in order to 
put it in a trellis oriented form, then the Wolff trellis is 
minimal. Indced, the trellis obtained in this way does 
not only minimize the number of vertices but also the 
state space profile. As for classical trellises, we define the 
state profile of the trellis of an (n, k) code by the (n + 1)- 
tuple (fo, £ii • • ■ 1 6») where & = log 2 |V;| where Vi is the 
i-th state space of the trellis. In other words, we are go- 
ing to prové that the Wolff trellis applied to a stabilizer 
set in trellis oriented form minimizes each £j individually. 
Without loss of generality, we now assume that the trellis 
is associated to the syndromc s = (0, ... 0). 

First, we define the trellis oriented form of a stabilizer 
set. For 1 < j < n — k, let c(j) and d(j) be respectively 
the position of the first (resp. last) component of Sj 
which is diffcrent from /. We say that the stabilizer set 
{Sj} is in trellis oriented form if and only if for all j (1) 
S c / ] = I for f > j + 1 and ^ Sf j) ; and (2) there 
is at most one j' ^ j such that b(j) = b(j') and in such 
case Sj ^ S^ ' . Note that any stabilizer set {Sj} can 
be put in trellis oriented form by indices permutations 
and group additions. 

Second, we show a lower bound on £j. For an (n, k) 
stabilizer code given by the stabilizer set {Sj}j, and for 
i G {0, . . . , n}, let C{ (the future subgroup) be the sub- 
group of S 1 - whose elements have their first i components 
equal to /; and let Cf (the past subgroup) be the sub- 
group of S^ whose elements have their last n — i compo- 
nents equal to /. We then have the following lemma. 

Lemma 1. £i>n + k- log 2 \Cf\ - log 2 |C/|. 

Proof. Let v G V;, and C v be the set of elements of S 1 - 
that correspond to a path in the trellis associated to the 
all-zero syndrome and passing through v. Let P and F 
be the set of all paths that go from vo to v (resp. from 
v to v„). Let Q be a fixed element of C v and qpqF its 
corresponding path in the trellis, where qp G P and qF £ 
F. By construction of the trellis, we have \C V \ = |P||F|. 
Note that any path p e P can be extended into a path 
pqF from Vq to v„. Since any path pqp can be associated 



to an element of Q + Cf, we have |P| < \Q + Cf\ = 
\Cf\. Similarly |F| < |C/|, which yiclds 2 n+k = ^ = 
Evev, |Cv| < \Vi\\Cf \\C{\. This gives the desired bound 
on by taking the logarithm. □ 

Finally we conclude with the following theorem. 

Theorem 1. The Wolff trellis achieves the previous 
bound on £j for each i when the stabilizer set is in trellis 
oriented form. 

Proof. Let 5 , s tart<í be the subset of operators of {Sj} 
which have at least one of their i-first components differ- 
ent from /. Let S^nd^i be the subset of operators of {Sj} 
which have all their last n — i components equal to /. It 
is straightforward to check that ^ = [S'stariK-tl ~ |Send<i|- 
Note that, Cf is the subgroup of G n with / on their 
n — i last components and orthogonal to all elements of 
S' s t a rt<i· This implies log 2 \Cf\ = 2i - |5 , s tart<i|· Using a 
similar argument, we get that log 2 \C{ \ = 2{n — i) — (n — 
k — |S'cnd<í|)· Adding these equalities, we conclude that 
6 + tog 2 ÏCf|+log 2 |C/|=n + fc. □ 

Using trellises of stabilizer codes. 

Min-Sum (Viterbi) algorithm. The Min-Sum algo- 
rithm is certainly one of the most widely employed al- 
gorithm that benefits from the trellis representation of 
classical codes. Here, we present a Min-Sum algorithm 
for stabilizer codes that computes the most likely error for 
memoryless Pauli channels given a measured syndrome s 
by using a trellis associated to the code. 

Consider an (n, k) stabilizer code with stabilizer 
set {Sj}j. Define the likclihood of P 6 G n as 
SlLi l°gP r ï(-P 1 )- Consider the n-section trellis for this 
quantum code associated to the syndrome s. The naming 
conventions for the vertices and edges are set as in pre- 
vious sections. For each edge e* of E^, define its weight 
wt(e l ) = — logPri(í(e ! )). By construction, the sum of 
weights along the path in the trellis that represents P is 
equal to the opposite of the likelihood. The task which 
consists in finding a most likely error E E G n with syn- 
drome s is thus equivalent to finding a lowest weight path 
(è 1 , . . . , è") in the trellis associated to s. 

This can be done by constructing recursively some sets 
Q of lowest weight error candidates. More precisely, Q 
contains couples (c, w) where c is a path issued from Vo 
that ends on a vèrtex of Vj and where w = wt(c). 

Min-Sum Algorithm 
Initialization: Co := {(vo,0)} and Q := for i < 1 
Main step: 
for i from 1 to n do 
for all v G Vi do 

Put in d the pair (c', wí(c')), where d is a path of 
minimum weight (ties are broken at random) among 
all paths that: (1) end in v; (2) have their i — 1 first 
vertices given by a path c of Q_i. 
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Note that the time complexity of this algorithm is linear 
in the number of vertices in the trcllis. 

Sum-Product algorithm. While the Min-Sum finds 
a most likely error compatible with the observed 
syndrome s, the Sum-Product aims at calculating 
marginal error probabilitics for physical qubits. That 
is, Pi(P) =Pr(error at qubit i = P\s) where P G G\. By 
dcfinition of the trellis, this probability is equal to the 
probability that a path (e l )^_ 1 from Vo to v„ is such that 
l{e l ) = P. 

The Sum-Product algorithm computes for each vèrtex 
v of the trellis associated with s a "forward" probability 
/(v) and a "backward" probability b(v). Both are then 
used to calculate the marginal probabilities Pi (P) . 

Sum-Product Algorithm 
Initialization: /(vq) := 1 and ò(v„) := 1; and /(v) := 



Above, (1) V~(v) is the set of vertices w in Vj_i which 
are adjacent to v; (2) V+(v) is the set of vertices w in 
V; + i which are adjacent to v; and (3) Ei(P) is the set of 
edges between Vj_i and Vi that bear the Pauli-label P. 

Once again, the practical relevance of this algorithm is 
due to the fact that its complexity is linear in the number 
of vertices in the trcllis. 



Computing the weight enumerator polynomial. The 
weight enumerator polynomial is a trivariate polynomial 
given by A(x, y, z) = J2i< u .v,w<n a u,v , w x u y v z w , where 
a u ,v,w = \{P S N(S) : \P\ X = u,\P\y = v,\P\ z = w}\ 
and where \P\e denotes the number of coordinates of 
P which are equal to E. It is possible to extract from 
A(x,y,z) a lot of useful information, e.g. bounds on the 
fidelity of the recovered state after decoding Q. As 
for classical codes, the weight enumerator can be com- 
puted with complexity linear in the number of vertices 
of the trellis. For this purpose intermediate polynomi- 
als A v (x,y, z) are calculated for each vèrtex of the trel- 
lis associated to s = (0, . . .0). Also define the polyno- 
mials Q I {x,y,z) = l, Q x (x,y,z)=x, Q Y (x,y,z)=ya,nd 
Qz(x,y,z) = z. 



Computation of A(x, y, z) 
Initialization: A VQ (x,y, z) = 1 
Main step: 
for i from 1 to n do 
for all v G Vi do 

A v (x,y,z) := E we v+(v) A ^(x,y, z)Qi (m) (x,y, z) 

A(x,y,z) := A Vn (x,y,z) 

Above, V~(v) is the set of vertices w in Vi_i which are 
adjacent to v. 

The proof of correetness for these three algorithms fol- 
lows from the correetness of the more general Min-Sum 
and Sum-Product algorithms presented in || . 

Conclusion. From a practical point of view, we have 
proposed a definition for the trellis of a stabilizer code 
together with two constructions and three algorithms 
that take advantage of this representation. Following 
the same path, several algorithms for classical codes run- 
ning on trellises can be generalized to the quantum case. 
Among them, estimation of noise parameters (or equal- 
ization) seems a promising avenue for enhancing the per- 
formance of quantum òptics fiber Communications using 
near-future quantum technology. 



and ò(v) := for all other vertices 
Forward pass : 
for i = 1 to 7i do 
for all v G V, do /(v) := £ we v r( v) /(w)P ri (/(wv)) 

for all v e Vi do /(v) := /(v)/F f 
Backward pass : 
for i = n — 1 down to do 

for all v G Vi do ò(v) = £ weV +( v ) &(w)Pr í+ i(Z(vw)) 

Bi := E veVl 6(v) 

for all v G V, do ò(v) := b{v)/Bi 
Final pass : 
for i = 1 to n do 

for all P G Gi do 

Pi(P) :=E V w e E ií P)/( v ) fc (w)Pr i (P) 
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